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1 Introduction 

In recent years the non-perturbative computation of the two point correlation func- 
tions of pure Yang-Mills theory have attracted a lot of attention. For what concerns 
the Landau gauge gluon propagator, the subject of this communication, lattice QCD 
simulations, Schwinger-Dyson equations and non-perturbative quantization of Yang- 
Mills theories provide essentially the same results. Recent reviews can be found in 
[T]. In this sense, we can claim to have a fair description of this two point function 
over the entire range of momentum. Here we review the large volume lattice simu- 
lations performed by one of the authors, its interpretation in terms of modeling the 
gluon propagator and we describe how the infrared propagator can be incorporated 
into an effective field theory model which approximates QCD at low energies. This 
effective theory connects gluon confinement with a gluon mass rUg with chiral sym- 
metry breaking. Further, the model predicts a particular simple relation between rUg 
and the light quark masses m^, i.e. that mg/nig is constant. This relation is tested 
using the solutions of the Schwinger-Dyson equations and found to be valid in the 
low energy regime below the 10% accuracy level. 

From the point of view of lattice simulations there are still some open questions. 
On the lattice, one numerically checks whether the selected Landau gauge configura- 
tions belong to the so-called Gribov region, which is the set of all transverse gauge 
connections with positive Faddeev-Popov operator. Within this set, one still finds 
Gribov copies, and it is not 100% well established if/how these additional copies 
influence the propagator in the (very deep) infrared region. However, previous simu- 
lations suggest that by choosing a different Gribov copy, the accompanying error lies 
typically within the statistical error of the propagator - see, for example, 
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Renormalized Gluon Propagator — p. = 3 GeV 
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Figure 1: DijP') for the largest volumes computed for each of the lattice spacings. 
The plot on the right shows D{p^) for the infrared region defined as j9 < 1 GeV. 



2 The Quenched Lattice Gluon Propagator 

In this section we report on the Landau gauge quenched lattice gluon propagator 

computed for large volumes La> 3.5 fm. The propagator has been computed previ- 
ously for huge lattice volumes, for the SU(3) gauge group up to -L a ^ 17 fm [1] and 
for SU(2) up to La ~ 27 fm, but using a ~ 0.2 fm [3]. The non-perturbative physics 
scale being ~ 1 fm, it is important to check these results by performing simulations 
at smaller lattice spacings and evaluate the corresponding finite volume effects. Here 
we report the computation of the propagator for a = 0.102 fm (/3 = 6.0), a = 0.0726 
fm (/3 = 6.2) and a = 0.0544 fm (/3 = 6.4) and for various volumes np to La ^ 8.2 
fm. The results reported in this section are preliminary. 

In Figured] we show the renormalized gluon propagator, computed with the differ- 
ent lattice spacings, for the largest physical volumes up to momenta p = 3 GeV and 
a zoom of the infrared region. For momenta above 1 GeV all data sets are, within 
errors, compatible. In the ultraviolet region, i.e. for p above ~ 2.8 GeV, all data sets 
are well described by the 1-loop inspired fit 

D{p') = , (2) 

where 7 = 13/22 is the gluon anomalous dimension. Indeed, for each set, the fits to ([2]) 



2 



Renormalized Gluoii Propagator - H = 3 GeV 
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Figure 2: Gluon mass from the infrared propagator (left) computed assuming a simple 
pole. The point on the right is the extrapolated mass to infinite volume. The right 
plot illustrates a typical fit. Note that the simple pole overestimates -D(O). 
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Table 1: -D(O) for the various simulations shown in Figure [H 



where used to renormalize the gluon propagator according to the MOM prescription 

Dr{p') = (3) 

The renormalization procedure is described in detail in [21 U] . Note, however, that in 
the present work, to renormalize, only the data for momenta above p ~ 2.8 GeV was 
included. In all the cases the "x^ jd.o.f. ~ 1. 

For the infrared region the data sets show finite volume and finite spacing effects. 
For example, the propagator computed using the simulation at /3 = 6.2, although 
having a smaller physical volume (5.84 fm)^, is below all the remaining data sets. For 
p = 0, the large statistical error hides the differences between the various -D(O). For 
the various simulations reported here, the corresponding -D(O) are given in Tabled] 
This numbers should be compared with the large volume,/? = 5.7, Wilson action SU(3) 
simulations performed by the Berlin-Moscow- Adelaide (BMA) group. Using the same 
definitions and setting the scale in the same way, i.e. from the string tension, the 
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BMA data reads D{0) = 8.68(37), 8.09(36), 7.59(56), 7.17(31) and 7.53(19) GeV'^ 
for La = 8.09, 11.76, 13.23, 14.70 and 16.17, respectively, given in fm. 

Qualitatively, the propagators computed with the different values of /3 are similar. 
The D{p^) from the /3 = 5.7 simulations are below all the remaining data sets. This 
suggest that, in order to quote a continuum propagator, one should extrapolate the 
lattice propagators to the infinite volume. Recall that, in principle, the renormal- 
ization procedure removes all the lattice spacing dependence. The extrapolation is 
better achieved if one uses a theoretical motivated functional form to describe -D(O) 
and extrapolates its parameters to the infinite volume - see also [6l [7] . 

Let us assume that the infrared propagator is described by a simple mass pole 

This functional form can only describe the propagator within a limited range of 
momenta. For example, a propagator described by does not violate positivity. 
Violation of positivity is a well established property of the non-perturbative gluon 
propagator. Anyway, at minimum, a fit to (jl]) defines an interval range where one 
can approximate -D(p^) by such a mass pole. The outcome of the fits as a function of 
the maximum range of momenta Pmax are 
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It follows that the lattice propagator can be described by a pole mass up to p ~ 500 
MeV with a gluon mass between 600 and 700 MeV. If one assumes that M{V) = 
Moo + Mi/L and fit the lattice data, then the infinite volume extrapolated mass is 
Moo = 634(40) MeV. The fit has a x^/d.o.f. = 1.4. The data for the different M and 
the extrapolation can be seen in Figure [2l Note that the simple pole (jl]) overestimates 
D{0), see right plot in Figure [21 and that is the reason why we do not provide the 
details of the extrapolation of Z and -D(O). 

The lattice data can be described by a propagator of type (jlj) if Z and M are 
functions of momentum, i.e. for 

A momentum dependent gluon mass together with a Z{p'^) were investigated in [8], 
where the same functional forms were used to fit the decoupling solutions of the 
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Figure 3: Gluon propagator and the fit to (j5]) using the definitions ([6]). Note that the 
functional form used here also overestimates -D(O). 
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Table 2: Fits of the lattice data to tree level RGZ propagator. La is given in fm and 
the mass scales and higher momenta in the fitting Pmax are given in power of GeV. 

Schwinger-Dyson equations. According to this work the lattice data can be described 

by 

Z{p')= ^ f'^, and M\p') = -^ (6) 

up to momenta Pmax = 4.2 GeV. The values of the fitted parameters for the largest 
physical volume are 

zo = 1.189(20), A = 1.842(39) GeV, r = 7.49(59) and mo = 671(9) MeV. 

The gluon data and the fit to ([5]) are reported in Figure [31 Note that, as in the case 
of simple pole mass, the fit overestimates -D(O). 
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Renormalized Gluon Propagator - |i = 3 GeV 
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Figure 4: Gluon propagator and the fits to ([7]) for tlie largest lattice volumes. Note 
that the RGZ reproduces well all the lattice data points, including -D(O). 



The refined Gribov-Zwanziger (RGZ) action is an improvement over the usual 
Faddeev-Popov quantization procedure for Yang-Mills theories, in the sense that it 
provides a better way to handle the problem of the Gribov copies by restricting 
the functional integration space to the so-called Gribov region. The RGZ action 
is renormalizable, in the perturbative sense, and introduces new auxiliary bosonic 
and fermionic fields. In what concerns the gluon propagator, the RGZ tree level 
propagator is given by 

2^ ^ p^ + M^ 

where is a mass scale related to the new auxiliary fields, m? is another mass scale 
related with the {A?) condensate and 7^ is the Gribov parameter. 7^ is not a free 
parameter but is fixed by the so-called horizon condition |9]. In the following we shall 
introduce the shorthand = 2g^N'~f'^ + M'^m'^ . The RGZ being a non-perturbative 
quantization for the Yang-Mills theories, one hopes that its tree level predictions 
provide a good description for the infrared. The propagator ([7]) can rewritten as 

D{p') = i-— with M^p") = + 2^7' AT 7^ 



In this sense, the RGZ action predicts a momentum dependent effective gluon mass 
which is essentially the functional form analyzed previously, i.e. M'^ijP') given by 
equation (|6]). The tree level expression for D{p'^) does not include the observed 
logarithmic corrections at high energies and, therefore, one expects ([7]) to deviate 
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Figure 5: Gluon propagator and the fit to ([7]) using the hnear extrapolated param- 
eters. The red fine was computed using the extrapolated {3 = 6.2 values, while the 
black line used the extrapolated f3 = 6.0 parameters. 



from the lattice data in the ultraviolet region. In order to extrapolate to the infinite 
volume and to have an estimate of the error on this extrapolation, we fit the following 
two sets of data to ([7]): (i) 48^ with La = 4.88 fm, 64^ with La = 6.50 fm and 80^ 
with La = 8.13 fm for /3 = 6.0; (ii) 48^ with La = 3.49 fm, 64^ with La = 4.65 fm and 
80^ with La = 5.81 fm for /3 = 6.2. The fits are summarized in Table |2l The lattice 
data and the fits for the largest physical volumes are reported in Figure |H Note that 
the RGZ propagator reproduces well all the lattice data, including -D(O). 

The infrared propagator can be extrapolated to the infinite volume if one assumes 
a linear dependence on l/{La). The extrapolations give 
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and the corresponding zero momentum value is 

r 6.90(93) GeV-2 from the /3 = 6.0 data set, 

^(0) = (9) 
[ 7.84(69) GeV"^ from the (5 = 6.2 data set. 

In the computation of -D(O), given the poor quality of the linear extrapolation for 
A^, we have used instead the fitted value from the largest physical volume. Figure [S] 
shows the lattice gluon propagator for the largest physical volumes computed using 
different lattice spacings and the extrapolated fits to ([7j) as described above. 
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In the RGZ propagator, the parameter rn? is related to the {A^). If one uses the 
figures from the extrapolations it follows 

2.71(41) GeV^ from the /3 = 6.0 data set, 

(/A2)3GeV = <| (10) 

3.21(48) GeV^ from the /3 = 6.2 data set. 



or 



2.45(38) GeV^ from the /3 = 6.0 data set, 

(/^')lOGeV= <! (11) 

2.90(43) GeV^ from the (3 = 6.2 data set. 
The values are shghtly below those reported in [5]. 

3 The Gluon Mass and Chiral Symmetry Breaking 

A gluon mass term in the QCD action is forbidden by gauge invariance and, therefore, 
a gluon mass rUg has to be generated dynamically. A non-vanishing mass means that 
the gluon field is short ranged. We should be careful not to attribute a physical mean- 
ing to this "massive" gluon, given the already mentioned positivity violation, which 
is a indication of the unphysical (confined) nature of the gluon. Besides providing 
the screening of the gluon, one may ask if there are additional implications of having 
nig 7^ 0. In this section, we show that, within an effective field theory for low energy 
QCD, a gluon mass is connected with chiral symmetry breaking, i.e. the theory either 
has rrig ^ and chiral symmetry is broken or chiral symmetry is restored and the 
gluon is a long range field. This section is based in the work [TO]. 

In QCD the fundamental fields are associated with quarks and gluons. However, 
to describe the low energy regime of QCD other fields can be included to define an 
effective theory. Let us assume that the non-perturbative physics is mainly associated 
with the gluon sector. Pure Yang-Mills theory has multi-gluon configurations as 
bound states. The simplest of these bound states is a two gluon state. Given that the 
gluon belongs to the adjoint representation, the two gluon state can be decomposed 
according to8(g)8 = 1©8©8©10©10©27. The lowest dimensional irrep is a singlet 
and can be identified with glueball states. The lightest glueball state has J^*-^ = O"*""*" 
and a predicted mass of ~ 1.7 GeV [llj. Such a mass scale is well above the usual 
low energy mass scales, ~ 1 GeV or lower, and therefore, from the point of view 
of an effective theory, one expects the singlet to play a minor role. The next lower 
dimensional representations are the two 8 representations. They distinguish amongst 
themselves because one of them is symmetric under interchange of the gluons, while 
the other one is antisymmetric. Of the 8 irreps only the symmetric representation 
can generate a scalar field, which can be written as 

r oc d^.^Fl^F'^^^ , (12) 
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where is the non-abehan Maxwell tensor. Of course, one can add to the above 
definition a quark contribution given by, for example, qt°-q, where t°- are the generators 
of the fundamental representation. Adding the two terms enables to estimate the 
contribution of quarks and gluons to the effective field, 

^-^^H^. (13) 

where A ~ ^qcd is a non-perturbative mass scale. Plugging into this the gluon 
condensate as{F^) = 0.04 GeV^ and the light quark condensate (qq) — (—270 MeV)^, 
it follows that the ratio gluon to quark content of 0" is around 7. 

Let us consider an effective theory which includes the gluon field ^4^, the quark 
fields qf , where / is a fiavor index, and an effective scalar field 0" that belongs to 
the adjoint representation of the SU(3) color group. In the following we will assume 
that the non-perturbative physics is contained in 0". Furthermore, being an effective 
field theory, it should describe hadronic physics only in the low energy regime and it 
does not need to be renormalizable. The effective Lagrangian reads 

^ = -Wf""" + j:%{^rD,-mf}qf 
^ f 



-G^^\q^tq 



- ^2 E ^"7^g] -FzY. ^"l^F)^^ r + h.c. (14) 

/ / 

where — d^j,-\-igT"'A^ is the covariant derivative, T'^ the SU(3) generators, uif the 
current quark mass associated with flavor /, Voct the effective potential associated 
with the scalar fleld. Cqf is the gauge flxing part of the Lagrangian and Cghost 
contains the ghost terms. The Lagrangian is gauge invariant, excepts for the Cgf 
term. The effective gauge coupling constant g parameterizes residual interactions and 
it should be a small number, i.e. one expects the theory can be treated perturbativcly. 
The new interactions with the scalar field, the terms proportional to G4, G5, Fi, F2 
and F3, where written assuming flavor independence of strong interactions. 

C includes the QCD Lagrangian and verifles the usual soft-pion theorems of chiral 
symmetry at low energy. The new interactions introduce new vertices, not present 
in the original QCD Lagrangian, which contribute to quark processes. Note that the 
only new quark color singlet operator mimics the ^Pq model describing OZI-allowed 
mesonic strong decays. 
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The 0" kinetic term couples to a quadratic gluon term through the operator 

]^g^ct)\T^T\,ct)<'AlA'^^. (15) 

If the scalar fields acquires a vacuum expectation value without breaking color sym- 
metry, i.e. 

(0") = and {(t)y^) = 5"^ (16) 
given that for the adjoint representation tr{T°'T^) = NcS""^, the gluon mass reads 

ml = N,g\\ (17) 

where = 3. From the definition it follows that i.e. v"^, and therefore 

the gluon mass is gauge invariant. The proof of gauge invariance follows from the 
transformation properties of 0°. 

In the same way, the operator [q q] (jy^-cj)"- shifts the quark masses giving rise to 
a constituent quark mass 

t\t2 _ -\ r< 

Mf = mf- {Nj - l)G,v' = - ^^^m]. (18) 

For light quarks, the constituent quark mass is given by the quark self energy which, 
in the model, is linked with the gluon mass. Note, for our definitions, that G5 is a 
negative number. If the constituent mass for the light quarks vanishes, chiral sym- 
metry should be broken dynamically, whereby the relation Mf oc in the effective 
model links chiral symmetry with a finite effective gluon mass. 

The quark condensate (g q) , an order parameter for chiral symmetry breaking, can 
be computed in the model as a function of the constituent quark mass, the gluon mass 
and the theory cut-off - see |T0] for details. Then, if one identifies the gluon mass with 
the mass measured from the lattice using a simple pole propagator, nig = 634 MeV, 
together with Mf = 330 MeV and {qq) = (—270 MeV)^, one is able to estimate some 
of the theory parameters: 

G 

CO = 879 MeV, gv = 366 MeV and -4 = -0.31 GeV"^ , 

9^ 

where uJ is the theory's cut-off. 

4 Testing a Model Prediction 

The effective model relates the constituent quark mass M and the gluon mass rrig 
through equation f|T8|) . For a vanishing current mass, equation f|T8l) predicts a constant 
value for the ratio M/m^, at least at tree level. This result can be tested looking 
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Figure 6: On the left hand side, the plot shows the quark masses from solving the 
fermionic SDE gap equation, using different ansatze for the quark-gluon vertex, and 
the squared gluon mass computed from quenched lattice simulations. Note that M(p^) 
depends slightly on the definition of the quark-gluon vertex. On the right hand side, 
the plot shows the ratio M/iv?. 



at the solutions of the Schwinger-Dyson equations. In the following we will use the 
results published in [I2]. For the gluon and ghost propagators, the authors used 
the results of lattice QCD simulations and solved the gap equation for a massless 
fermion. The calculation does not take into account fermion loops and can be viewed 
as a quenched approximation. 

In [12], the fermionic gap equation was solved for two different ansatze for the 
quark-gluon vertex, a non-Abelian improved version of the Ball-Chiu vertex and an 
improved version of the Curtis-Pennington vertex. The choice of vertex leads to 
slightly different quark mass. In order to distinguish, the results of the Ball-Chiu 
vertex will be referred as BC, while the results from using the Curtis-Pennington 
vertex will be referred as CP. Figure |6] shows M computed from the Schwinger-Dyson 
equations for the different vertex ansatze, together with function of and, on 

the right hand side, the ratio M/m'^. The plots shows that M/m'^ increases slightly. 
If one looks at the maximal momentum range where the lattice gluon propagator can 
be fitted by a simple pole, i.e. if one compares the ratios up to momenta p ~ 0.5 
GeV, then M/m?g changes by less than 8%, relative to its zero momentum value, when 
using the BC quark-gluon vertex and less than 10% when using the CP vertex. 

5 Results and Conclusions 

We have currently a fair description of the gluon propagator over all momentum 
ranges. To extract the various parameters modeling the propagator, it would be 
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desirable to perform a high statistic and large volume simulation. 

The results of lattice simulations and Schwinger-Dyson equations show that the 
gluon propagator behaves as a dynamically massive gauge boson in the infrared region, 
see also the discussion in [T3] and references therein, and, the effective model sketched 
here, shows a connection between the gluon mass and chiral symmetry breaking. 
Comparing the tree level mass ratio prediction with the solutions of the Schwinger- 
Dyson equations we found good agreement in the low energy regime. 
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